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We present a measurement of the inclusive electron spectrum in B — ► X u eu decays near the 
kinematic limit for B — ► X a ev transitions, using a sample of 88 million BB pairs recorded by the 
BABAR detector at the T(4S) resonance. Partial branching fraction measurements are performed 
in five overlapping intervals of the electron momentum; for the interval of 2.0 - 2.6 GeV/ c we obtain 
AB(B -> X u ev) = (0.572 ± 0.041 s t a t ± 0.065 SHSt ) x 10~ 3 . Combining this result with shape function 
parameters extracted from BABAR measurements of moments of the inclusive photon spectrum in 
B — > X s 7 decays and moments of the hadron mass and lepton energy spectra in B — > X c lv decays 
we determine \V u b\ = (4.44 ± 0.25 M p ^o.38SF ^ 0.22 t heor V ) x 10~ 3 . Here the first error represents 
the combined statistical and systematic experimental uncertainties of the partial branching fraction 
measurement, the second error refers to the uncertainty of the determination of the shape function 
parameters, and the third error is due to theoretical uncertainties in the QCD calculations. 

PACS numbers: 13.20.He, 12.15.Hh, 12.38.Qk, 14.40.Nd 



I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ment V u b , the coupling of the b quark to the u quark, is 
a fundamental parameter of the Standard Model. It is 
one of the smallest and least known elements of the CKM 
matrix. With the increasingly precise measurements of 
decay-time-dependent CP asymmetries in .B-meson de- 
cays, in particular the angle (3 [1, 2], improved measure- 
ments of the magnitude of V u t, will allow for stringent 
experimental tests of the Standard Model mechanism for 
CP violation [3]. This is best illustrated in terms of 
the unitarity triangle, the graphical representation of the 
unitarity condition for the CKM matrix, for which the 
length of the side that is opposite to the angle /? is pro- 
portional to \V u b\- 

The extraction of \V u b\ is a challenge, both theoreti- 
cally and experimentally. Experimentally, the principal 
challenge is to separate the signal B —* X u ev decays from 
the 50 times larger B — > X c ev background. This can 
be achieved by selecting regions of phase space in which 
this background is highly suppressed. In the rest frame 
of the B meson, the kinematic endpoint of the electron 
spectrum is ~ 2.3GcV/c for the dominant B — > X c ev 
decays and <~ 2.6GeV/c for B — > X u ev decays. Thus 
the spectrum above 2.3 GeV/c is dominated by electrons 
from B — > X u ev transitions. This allows for a relatively 
precise measurement, largely free from BB background, 
in a 300 MeV/c interval that covers approximately 10% 
of the total electron spectrum for charmless semileptonic 
B decays. In the T(45) rest frame, the finite momenta of 
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the B mesons cause additional spread of the electron mo- 
menta of ~ 200 MeV/c, extending the endpoints to higher 
momenta. 

The weak decay rate for b —* uev can be calculated 
at the parton level. It is proportional to \V u b\ 2 and 
ml, where nib refers to the 6-quark mass. To relate 
the semileptonic decay rate of the B meson to \V u b\, 
the parton-level calculations have to be corrected for 
perturbative and non-perturbative QCD effects. These 
corrections can be calculated using various techniques: 
heavy quark expansions (HQE) [4] and QCD factoriza- 
tion [5]. Both approaches separate perturbative from 
non-perturbative expressions and sort terms in powers 
of 1/mb. HQE is appropriate for the calculations of to- 
tal inclusive B decay rates and for partial B decay rates 
integrated over sufficiently large regions of phase space 
where the mass and momentum of the final state hadron 
are large compared to Aqcd- QCD factorization is bet- 
ter suited for calculations of partial rates and spectra 
near the kinematic boundaries where the hadronic mass 
is small. In this region the spectra are affected by the 
distribution of the 6-quark momentum inside the B me- 
son [5], which can be described by a structure or shape 
function (SF), in addition to weak annihilation and other 
non-perturbative effects. Extrapolation from the limited 
momentum range near the endpoint to the full spectrum 
is a difficult task, because the SF cannot be calculated. 
To leading order, the SF should be universal for all b — > q 
transitions (here q represents a light quark) [6, 7]. Sev- 
eral functional forms for the SF, which generally depend 
on two parameters related to the mass and kinetic en- 
ergy of the 6-quark, A or mj, and Ai or u, 2 , have been 
proposed. The values and precise definitions of these pa- 
rameters depend on the specific ansatz for the SF, the 
mass renormalization scheme, and the renormalization 
scale chosen. 

In this paper, we present a measurement of the inclu- 
sive electron momentum spectrum in charmless semilep- 
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tonic B decays, averaged over charged and neutral B 
mesons, near the kinematic endpoint. We report mea- 
surements of the partial branching fractions in five over- 
lapping momentum intervals. The upper limit is fixed at 
2.6 GeV/c, while the lower limit varies from 2.0 GeV/c 
to 2.4 GeV/c. By extending the interval for the signal 
extraction down to 2.0 GeV/c, we capture about 25% of 
the total signal electron spectrum, but also much larger 
B — > X c ev backgrounds. Inclusive measurements of \V u b\ 
have been performed by several experiments operating at 
the T(4S*) resonance, namely ARGUS [8], CLEO [9, 10], 
BABAR [11], and Belle [12], and experiments operating 
at the Z° resonance, namely L3 [13], ALEPH [14], DEL- 
PHI [15], and OPAL [16]. This analysis is based on a 
method similar to the one used in previous measurements 
of the lepton spectrum near the kinematic endpoint [8, 9]. 
The results presented here supersede those of the prelim- 
inary analysis reported by the BABAR Collaboration [11]. 

The extraction of \V u b\ relies on two different the- 
oretical calculations of the differential decay rates for 
B — ► X u ev and B — > X s -y: the original work by DeFazio 
and Neubert (DN) [17], and Kagan and Ncubert [18], and 
the more comprehensive recent calculations by Bosch, 
Lange, Neubert, and Paz (BLNP) [19-24]. 

The DN calculations allow for the extrapolation of the 
observed partial B — > X u ev decay rate above a certain 
electron momentum to the total inclusive B —>■ X u ev 
decay rate using the measured SF parameters and a sub- 
sequent translation of the total decay rate to \V u b\. The 
theoretical uncertainties on the rate predictions are esti- 
mated to be of order 10-20%. 

The BLNP authors have presented a systematic treat- 
ment of the SF effects, incorporated all known corrections 
to the differential decay rates, and provided an interpo- 
lation between the HQE and the SF regions. They have 
also performed a detailed analysis of the theoretical un- 
certainties. The calculations directly relate the partial 
decay rate to \V u b\. While the calculations by BLNP 
are to supersede the earlier work by DN, we use both 
approaches to allow for a direct comparison of the two 
calculations, and also a comparison with previous mea- 
surements based on the DN calculations. We adopt the 
SF parameters extracted by the BABAR Collaboration: 
for the DN method we rely on the photon spectrum in 
B — ► X s j decays [25]; for the more recent BLNP method, 
we also use SF parameters derived from the photon spec- 
trum, its moments, the hadron-mass and lepton-energy 
moments in inclusive B — ► X c lv decays [26], and the 
combination of all moments measured by the BABAR Col- 
laboration [27]. 



II. DATA SAMPLE, DETECTOR, AND 
SIMULATION 

The data used in this analysis were recorded with the 
BABAR detector at the PEP-II energy-asymmetric e + e~ 
collider. The data sample of 88 million BB events, cor- 



responding to an integrated luminosity of 80.4 fb~ , was 
collected at the T(45) resonance. An additional sam- 
ple of 9.5 fb _1 was recorded at a center-of-mass (cm.) 
energy 40 McV below the T(45) resonance, i.e. just be- 
low the threshold for BB production. This off-resonance 
data sample is used to subtract the non-BB contribu- 
tions from the data collected on the T(45) resonance. 
The relative normalization of the two data samples has 
been derived from luminosity measurements, which are 
based on the number of detected \i + [i~ pairs and the 
QED cross section for e+e~ — ► production, ad- 

justed for the small difference in center-of-mass energy. 

The BABAR detector has been described in detail else- 
where [28]. The most important components for this 
study are the charged-particle tracking system, consist- 
ing of a five-layer silicon detector and a 40-layer drift 
chamber, and the electromagnetic calorimeter assembled 
from 6580 CsI(Tl) crystals. These detector components 
operate in a 1.5-T solenoidal magnetic field. Electron 
candidates are selected on the basis of the ratio of the 
energy detected in the calorimeter to the track momen- 
tum, the calorimeter shower shape, the energy loss in the 
drift chamber, and the angle of the photons reconstructed 
in a ring-imaging Chcrcnkov detector. 

The electron identification efficiency and the probabil- 
ities to misidentify a pion, kaon, or proton as an electron 
have been measured [29] as a function of the laboratory 
momentum and angles with clean samples of tracks se- 
lected from data. Within the acceptance of the calorime- 
ter, defined by the polar angle in the laboratory frame, 
—0.72 < cos 9i a b < 0.92, the average electron identifica- 
tion efficiency is 92%. The average hadron misidentifica- 
tion rate is about 0.1%. 

We use Monte Carlo (MC) techniques to simulate the 
production and decay of B mesons, and the detector 
response [30], to estimate signal and background effi- 
ciencies, and to extract the observed signal and back- 
ground distributions. The simulated sample of generic 
BB events exceeds the BB data sample by about a fac- 
tor of three. 

Information from studies of selected control data sam- 
ples on efficiencies and resolutions is used to improve the 
accuracy of the simulation. Comparisons of data with 
the MC simulations have revealed small differences in 
the tracking efficiencies, which have been corrected for. 
No significant impact of non-Gaussian resolution tails has 
been found for high momentum tracks in the endpoint re- 
gion. The MC simulations include radiative effects such 
as bremsstrahlung in the detector material and QED ini- 
tial and final state radiation [31]. Adjustments for small 
variations of the beam energy over time have also been 
included. 

In the MC simulations the branching fractions for 
hadronic B and D decays are based on values reported 
in the Review of Particle Physics [32]. The simulation 
of charmless semileptonic decays, B — > X u ev, is based 
on a heavy quark expansion to 0(a s ) [17]. This calcula- 
tion produces a continuous spectrum of hadronic states. 
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Electron Momentum (GeV/c) 

FIG. 1: MC-generated electron momentum spectra for vari- 
ous charmless semileptonic B decays: B — > neis, B — > pev, 
B — > cuev, B — > ryej/, B — > the sum of B-meson de- 

cay modes to non-resonant and higher-mass resonance states 
{X^ev), and the sum of all decay modes (All). The spectra 
are normalized to a total rate of 1.0. 



The hadronization of X u with masses above 2m n is per- 
formed by JETSET [33]. The motion of the b quark 
inside the B meson is implemented with the SF parame- 
terization given in [17]. Three-body decays to low-mass 
hadrons, (X u — n, p, ui, r), rf), are simulated separately 
using the ISGW2 model [34] and mixed with decays to 
non-resonant and higher mass resonant states X* 7 so that 
the cumulative distributions of the hadron mass, the mo- 
mentum transfer squared, and the electron momentum 
reproduce the HQE calculation as closely as possible. 
The generated electron spectrum is reweighted to accom- 
modate variations due to specific choices of the SF pa- 
rameters. 

The MC-generated electron-momentum distributions 
for B — ► X u ev decays are shown in Fig. 1, for individual 
decay modes and for their sum. Here and throughout the 
paper, the electron momentum and all other kinematic 
variables are measured in the T(AS) rest frame, unless 
stated otherwise. Above 2 GeV/ c, the principal signal 
contributions are from decays involving the light mesons 
7r, p, and u>, and also some higher mass resonant and non- 
resonant states X*. 

For the simulation of the dominant B — > X c ev decays, 
we have chosen a variety of models. For B — > Dev and 
B — > D*ev decays we use parameterizations [35-37] of 
the form factors, based on heavy quark effective theory 
(HQET). Decays to pseudoscalar mesons are described by 
a single form factor Fd (w)/ Fd (1) = 1 — p D (w — 1) , where 
the variable w is the scalar product of the B and D meson 
four- vector velocities and is equal to the relativistic boost 
of the D meson in the B meson rest frame. The linear 
slope p 2 D has been measured by the CLEO [38] and Belle 
[39] Collaborations. We use the average value, p 2 D — 



0.72±0.12. The differential decay rate for B — > D*ev can 
be described by three amplitudes, which depend on three 
parameters: p 2 , R\, and R2. We adopt values recently 
measured by BABAR [40]: p 2 = 0.769 ± 0.043 ± 0.032, 
Ri = 1.328 ± 0.060 ± 0.025, and R 2 = 0.920 ± 0.048 ± 
0.013. Here the parameter p 2 is the slope assuming a 
linear dependence of the form factor on the variable w. 
The quoted errors reflect the statistical and systematic 
uncertainties. 

We use the ISGW2 [34] model for various decays to 
higher-mass D** resonances. We have adopted a pre- 
scription by Goity and Roberts [41] for the non-resonant 
B — > D^irev decays. 

The shapes of the MC-generated electron spectra for 
individual B — ^ X c ev decays are shown in Fig. 2. Above 
2 GeV/c the principal background contributions are from 
decays involving the lower-mass charm mesons, D* and 
D. Higher-mass and non-resonant charm states are ex- 
pected to contribute at lower electron momenta. The 
relative contributions of the individual B — > X c ev decay 
modes are adjusted to match the data by a fit to the 
observed spectrum (see below). 




0.0 1.0 2.0 3.0 

Electron Momentum (GeV/c) 



FIG. 2: MC-generated electron momentum spectra for various 
B —* X c eu decay modes: B — > Dev, B — > D*eu, B —* D**ev, 
B — > D^Tvev, and B — > X u eu, and the sum of all decay 
modes (All). The signal B — > X u eu spectrum is shown for 
comparison. The spectra are normalized to a total rate of 
1.0. 



III. ANALYSIS 
A. Event Selection 

We select events with a semileptonic B decay by re- 
quiring an electron with momentum p e > 1.1 GeV/c. To 
reject electrons from the decay J/ip — > e + e~, we com- 
bine the electron candidate with any second electron of 
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opposite charge and reject the combination, if the invari- 
ant mass of the pair falls in the interval 3.00 < m ee < 
3.15GeV/c 2 . 

To suppress background from non-BB events, primar- 
ily low-multiplicity QED (including t + t~ pairs) and 
e+e~ — ► qq processes (here q represents any of the u, d, s 
or c quarks), we veto events with fewer than four charged 
tracks. We also require that the ratio of the second to the 
zeroth Fox- Wolfram moment [42], IZ2, not exceed 0.5. 7?-2 
is calculated including all detected charged particles and 
photons. For events with an electron in the momentum 
interval of 2.0 to 2.6 GeV/c, these two criteria reduce the 
non-BB background by a factor of about 6, while the loss 
of signal events is less than 20%. 

In semileptonic B decays, the neutrino carries siz- 
able energy. In events in which the only unde- 
tected particle is this neutrino, the neutrino four- 
momentum can be inferred from the missing momen- 
tum, p m i SS = (E m i SS ,p miss ), the difference between the 
four-momentum of the two colliding-beam particles, and 
the sum of the four- momenta of all detected particles, 
charged and neutral. To improve the reconstruction of 
the missing momentum, we impose a number of require- 
ments on the charged and neutral particles. Charged 
tracks are required to have a minimum transverse mo- 
mentum of 0.2 GcV/c and a maximum momentum of 
10 GcV/c in the laboratory frame. Charged tracks are 
also restricted in polar angle to —0.82 < cos 9i a b < 0.92 
and they are required to originate close to the beam- 
beam interaction point. The individual photon energy 
in the laboratory frame is required to exceed 30 MeV. 
The selection of semileptonic B decays is enhanced by 
requiring |p m j ss | > 0.5GeV/c, and that p m i SS points into 
the detector fiducial volume, \ cos 8 m i SS \ < 0.9, thereby 
effectively reducing the impact of particle losses close to 
the beams. Furthermore, since in semileptonic B decays 
with a high-momentum electron, the neutrino and the 
electron are emitted preferentially in opposite directions, 
we require that the angle Aa between these two particles 
fulfill the condition cosAa < 0.4. These requirements 
for the missing momentum reduce the continuum back- 
ground from QED processes and e + e~ — » qq production 
by an additional factor of 3, while the signal loss is less 
than 20%. 

The stated selection criteria result in an efficiency (in- 
cluding effects of bremsstrahlung) of 35 — 50% for select- 
ing B — > X u ev decays; its dependence on the electron 
momentum is shown in Fig. 3. 



B. Background Subtraction 

The spectrum of the highest momentum electron in 
events selected by the criteria described above is shown 
in Fig. 4a, separately for data recorded on and below the 
T(4S) resonance. The data collected on the T(45) reso- 
nance include contributions from BB events and non-Si? 
background. The latter is measured using off-resonance 
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FIG. 3: Selection efficiency for events with B — > X u ev de- 
cays as a function of the electron momentum. The error bars 
represent the statistical errors. 



data, collected below BB production threshold, and us- 
ing on-resonance data above 2.8 GeV/c, i.e., above the 
endpoint for electrons from B decays. The BB back- 
ground to the B — ► X u ev spectrum is estimated from 
MC simulation, with the normalization of the individual 
contributions determined by a fit to the total observed 
spectrum. 



1. Non-BB Background 

To determine the non-£?_B background we perform a 
X 2 fit to the off-resonance data in the momentum inter- 
val of 1.1 to 3.5 GeV/c and to on-resonance data in the 
momentum interval of 2.8 to 3.5 GeV/c. Since the cm. 
energy for the off-resonance data is 0.4% lower than for 
the on-resonance data, we scale the electron momenta for 
the off-resonance data by the ratio of the cm. energies. 

The relative normalization for the two data sets is 



tl = 



SOFF jLoNdt 

sqn jLoFFdt 



= 8.433 ±0.004 ±0.021, 



where s and fLdt refer to the cm. energy squared and 
integrated luminosity of the two data sets. The statis- 
tical uncertainty of tl is determined by the number of 
detected pairs used for the measurement of the in- 

tegrated luminosity; the systematic error of the ratio is 
estimated to be 0.25%. 

The x 2 f° r the fit to the non-BB events is defined as 
follows, 



= E 



(f(a,Pi) ~ r L n t ) 2 



if(a,Pj) 



N 3 ? 



(1) 
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FIG. 4: (color online) Electron momentum spectra in the 
T(45) rest frame: (a) on-resonance data (open circles - blue), 
scaled off-resonance data (solid circles - green); the solid line 
shows the result of the fit to the non-BB events using both on- 
and off-resonance data; (b) on-resonance data after subtrac- 
tion of the fitted non-BB background (triangles - blue) com- 
pared to simulated BB background that is adjusted by the 
combined fit to the on- and off-resonance data (histogram); 
(c) on-resonance data after subtraction of all backgrounds 
(linear vertical scale, data points - red), compared to the sim- 
ulated B — > X u eu signal spectrum (histogram); the error bars 
indicate errors from the fit, which include the uncertainties in 
the fitted scale factors for non-SB and X c ev backgrounds. 
The shaded area indicates the momentum interval for which 
the on-resonance data are combined into a single bin for the 
purpose of reducing the sensitivity of the fit to the shape of 
the signal spectrum in this region. 



Here rii and Nj refer to the number of selected events 
in the off- and on-resonance samples, for the i-th or j-th. 
momentum bin (pj > 2.8GeV/c), and a is the set of free 
parameters of the fit. For the function approximating 
the momentum spectrum, we have chosen an exponential 
expression of the form 



/(a,p) = ai + exp(a 2 + a 3 p + cup 2 ) 



(2) 



The fit describes the data well: \ 2 = 70 for 58 degrees of 
freedom. Above 2.8 GeV/c, we observe (36.7 ±0.2) x 10 3 
events in the on-resonance data, compared to the fitted 
number of (36.6 ± 0.2) x 10 3 events. 



2. BB Background 

The electron spectrum from B-meson decays is com- 
posed of several contributions, dominated by the vari- 



ous semileptonic decays. Hadronic B decays contribute 
mostly via hadron misidentification and secondary elec- 
trons from decays of D, J/tp, and ip(2S) mesons. 

We estimate the total background by fitting the ob- 
served inclusive electron spectrum to the sum of the sig- 
nal and individual background contributions. For the 
individual signal and BB background contributions, we 
use the MC simulated spectra, and treat their relative 
normalization factors as free parameters in the fit. The 
non-BB background is parameterized by the exponential 
functions f(a,Pi), as described above. We expand the x 2 
definition as follows, 



X 



rLni) 2 



E 



(f(a,Pj) + S(b,pj) 



N, 



a 



(3) 



jMC 



where the first sum is for the off-resonance data and the 
second sum for the on-resonance data. The BB electron 
spectrum is approximated as S(b,pj) = J2k^k9k(Pj), 
where the free parameters bk are the correction factors 
to the MC default branching fractions for the six indi- 
vidual contributions gk{Pj) representing the signal B — » 
X u ev decays, the background B — ► Deis, B — > D*ev, 
B -» D**ev, B -> DW ireu decays, and the sum of other- 
background events with electrons from secondary decays 
or misidentified hadrons. dj mc is the statistical error of 
the number of simulated events in the j-th bin. The mo- 
mentum spectra gk(Pj) are histograms taken from MC 
simulations. 



3. Fit to Inclusive Spectra 

The fit is performed simultaneously to the on- and off- 
resonance electron momentum spectra in the range from 
1.1 to 3.5 GeV/c, in bins of 50MeV/c. The lower part of 
the spectrum determines the relative normalization of the 
various background contributions, allowing for an extrap- 
olation into the endpoint region above 2.0 GeV/c. To re- 
duce a potential systematic bias from the assumed shape 
of the signal spectrum, we combine the on-resonance data 
for the interval from 2.1 to 2.8 GeV/c into a single bin. 
The lower limit of this bin is chosen so as to retain the 
sensitivity to the steeply falling BB background distri- 
butions, while containing a large fraction of the signal 
events in a region where the background is low. The fit 
results are insensitive to changes in this lower limit in the 
range of 2.0 to 2.2 GeV/c. The number of signal events 
in a given momentum interval is taken as the excess of 
events above the fitted background. 

The observed spectra, the fitted non-BB and BB back- 
grounds and the signal are shown and compared to MC 
simulations in Fig. 4. The fit has a % 2 of 96 for 73 de- 
grees of freedom. Above 2.3 GeV/c, the non-BB back- 
ground is dominant, while at low momenta the semilep- 
tonic BB background dominates. Contributions from 
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TABLE I: Summary of the signal extraction: the number of events (in units of 10 3 ) for the total sample, the principal background 
contributions, and the remaining signal, as well as the signal efficiencies, for five intervals of the electron momentum. The 
errors listed are statistical, including the uncertainties in the fitted scale factors for the non-BB and X c ev backgrounds. The 
error values of 0.00 represent errors of less than 0.005. 



Ap ( GeV/c) 


2.0 - 2.6 


2.1 - 2.6 


2.2 - 2.6 


2.3 - 2.6 


2.4 - 2.6 


Total sample 


609.81 ± 0.78 


295.76 ± 0.54 


133.59 ± 0.37 


65.48 ± 0.26 


35.38 ± 0.19 


ly\Jll-U>U UdLKgl ULII1U 


1 49 ^8 4- fl f\\ 


i 90 + n 4s 


74 Sfi + n "\(\ 
li.ou zn u.ou 


uu.-LO zn \j.z,o 


9Q Qfi 4- fl 1 (S 
iy.yu ziz u.iu 


X c ev background 


416.22 ± 2.52 


157.17 ± 1.29 


38.82 ± 0.47 


4.00 ± 0.10 


0.09 ± 0.01 


J/tp and %1>{2S) 


6.17 ± 0.14 


4.00 ± 0.10 


2.33 ± 0.06 


1.17 ± 0.04 


0.47 ± 0.02 


Other e background 


1.61 ± 0.05 


0.62 ± 0.02 


0.24 ± 0.01 


0.08 ± 0.01 


0.03 ± 0.00 


7r mis-identification 


1.34 ± 0.04 


0.98 ± 0.03 


0.64 ± 0.02 


0.34 ± 0.02 


0.10 ± 0.01 


K mis-identification 


0.47 ± 0.02 


0.26 ± 0.01 


0.13 ± 0.01 


0.05 ± 0.01 


0.01 ± 0.00 


Other mis-identification 


0.27 ± 0.01 


0.15 ± 0.01 


0.08 ± 0.01 


0.04 ± 0.01 


0.02 ± 0.00 


X u ev background 


1.62 ± 0.10 


0.66 ± 0.05 


0.20 ± 0.02 


0.03 ± 0.01 


0.01 ± 0.00 


X u ev signal 


39.72 ± 2.70 


26.72 ± 1.49 


16.31 ± 0.71 


9.64 ± 0.38 


4.70 ± 0.25 


X u ev efficiency (%) 


42.1 ± 0.3 


41.2 ± 0.4 


40.2 ± 0.5 


39.5 ± 0.7 


37.9 ± 1.0 



hadron misidcntification arc small, varying from 6% to 
4% as the electron momentum increases. The theoreti- 
cal prediction for the signal B —> X u ev spectrum based 
on the BLNP calculations uses SF parameters extracted 
from the combined fit [27] to the moments measured by 
the BABAR Collaboration. 

The fitting procedure was chosen in recognition of the 
fact that currently the branching fractions for the indi- 
vidual B — > X c iv decays are not well enough measured to 
perform an adequate background subtraction. The MC 
simulation takes into account the form factor and angu- 
lar distributions for the B — > Dev and B — ► D*ev de- 
cays. For decays to higher-mass mesons, this information 
is not available. As a result, wc do not consider this fit 
as a viable method of measuring these individual branch- 
ing fractions. Nevertheless, the fitted branching fractions 
agree reasonably well with the measured branching frac- 
tions [32]. For the decays to higher- mass states, the abil- 
ity of the fit to distinguish between decays to D**ev and 
D*irev is limited. The sum of the two contributions, 
however, agrees with current measurements [32]. 

Tabic I shows a summary of the data, principal back- 
grounds and the resulting signal. The errors are statisti- 
cal, but for the non-SB and X c ev background they in- 
clude the uncertainties of the fitted parameters. The data 
are shown for five overlapping signal regions, ranging in 
width from 600 to 200 MeV/c. We choose 2.6 GcV/c as 
the common upper limit of the signal regions because at 
higher momenta the signal contributions are very small 
compared to the non-BB background. As the lower limit 
is extended to 2.0 GcV/c, the error on the BB back- 
ground subtraction increases. 



IV. SYSTEMATIC ERRORS 

A summary of the systematic errors is given in Table II 
for five intervals in the electron momentum. The prin- 
cipal systematic errors originate from the event selection 



and the background subtraction. The uncertainty in the 
event simulation and its impact on the momentum de- 
pendence of the efficiencies for signal and background are 
the experimental limitations of the current analysis. The 
second largest source of uncertainties is the estimate of 
the BB background derived from the fit to the observed 
electron spectrum, primarily due to the uncertainties in 
the simulated momentum spectra of the various contribu- 
tions. In addition, there are relatively small corrections 
to the momentum spectra due to variations in the beam 
energies, and radiative effects. 



A. Detection and Simulation of B — » X u ev Decays 

The selection efficiency for B — > X u ev decays is deter- 
mined by MC simulation. We include in the uncertainty 
of the signal spectrum not only the uncertainty in the 
simulation of the detector response, but also the uncer- 
tainty in the simulation of the momentum and angular 
distributions of the electron, as well as the hadrons and 
the neutrino. 



1. Detector related uncertainties 

For a specific model of the signal decays there are three 
major factors that determine the efficiency: the track re- 
construction for the electron, the electron identification, 
and losses due to the detector acceptance and event se- 
lection. 

The uncertainty in the tracking efficiency has been 
studied in detail and is estimated to be ~ 0.7% per 
track. The average identification efficiency for electrons 
with momentum above 1.0 GeV/c is estimated to be on 
average 92% [29], based on large control samples of ra- 
diative Bhabha events and two-photon interactions. In 
BB events the actual efficiencies are slightly lower due 
to higher track and photon multiplicity. This difference 
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TABLE II: Summary of the relative systematic errors (%) on 
the partial branching fraction measurements for B — > X u ev 
decays, as a function of p mm , the lower limit of the signal 
momentum range. The common upper limit is 2.6 GeV/c. 
The sensitivity of the signal extraction to the uncertainties in 
the SF parameters is listed as an additional systematic error, 
separately for the four sets of SF parameters. 



p min ( GeV/ c) 


2.0 


2.1 


2.2 


2.3 


2.4 


Track finding efficiency 


0.7 


0.7 


0.7 


0.7 


0.7 


Electron identification 


1.4 


1.4 


1.3 


0.9 


0.8 


Event selection efficiency 


6.8 


6.7 


6.1 


5.5 


7.9 


Non-SB background 


2.4 


2.5 


2.4 


2.5 


2.3 


J/ip and tp(2S) background 


0.9 


0.8 


0.8 


0.6 


0.5 


B — > D*lv form factor 


2.4 


2.3 


2.0 


1.3 


0.5 


B — ► Dlv form factor 


0.7 


0.9 


0.8 


0.2 


0.4 


B — > D**cis spectrum 


2.8 


2.5 


2.4 


0.9 


0.7 


Other e background 


0.5 


0.3 


0.2 


0.1 


0.1 


B — > X u ev background 


1.1 


0.6 


0.3 


0.1 


0.0 


7r mis-identification background 


0.8 


0.9 


0.9 


0.8 


0.5 


K mis— identification background 


0.4 


0.4 


0.3 


0.2 


0.1 


Other hadron mis-identification 


0.2 


0.2 


0.2 


0.1 


0.1 


B movement 


1.3 


1.7 


1.5 


0.6 


0.1 


Bremsstrahlung and FSR 


1.0 


1.2 


1.2 


0.9 


0.9 


N B g normalization 


1.1 


1.1 


1.1 


1.1 


1.1 


Total experimental error 


8.8 


8.6 


7.9 


6.6 


8.5 


B — > X u eu spectrum 












Xs'y SF, fit to spectrum 


6.0 


3.5 


1.6 


0.3 


0.1 


X s 7 SF, fit to moments 


11.3 


6.7 


3.1 


0.6 


0.1 


X c ev SF, fit to moments 


13.3 


8.6 


4.0 


0.8 


0.0 


SF, combined fit to moments 


7.2 


4.8 


2.3 


0.5 


0.0 


Total systematic error 












Xsj SF, fit to spectrum 


10.7 


9.3 


8.1 


6.6 


8.5 


Xs'j SF, fit to moments 


14.3 


10.9 


8.5 


6.6 


8.5 


X c ev SF, fit to moments 


15.9 


12.2 


8.9 


6.6 


8.5 


SF, combined fit to moments 


11.4 


9.8 


8.2 


6.6 


8.5 



decreases gradually from about 2.5% at 1.0 GeV/c to less 
than 0.8% at 2.0 GeV/c and above. We add in quadra- 
ture 50% of this observed difference to the statistical and 
systematic errors from the control samples. We assess 
the impact of this momentum-dependent uncertainty on 
the observed electron spectrum for both signal and back- 
ground (see below). 



2. Uncertainties in the signal spectrum 

The momentum distribution of the signal electrons is 
not precisely known because many of the exclusive de- 
cay modes that make up the total inclusive B — ► X u ev 
decays are still unobserved or poorly measured due to 
small event samples, and the form factors for most of 
the observed exclusive decay modes are not measured. 
To evaluate the sensitivity of the signal efficiency to the 
decay multiplicity and the shape of the momentum spec- 
trum, we independently vary the relative contributions of 
the different decay modes by their current experimental 
uncertainties. We observe changes in the signal yield of 



less than 3.0% for the spectrum above 1.1 GeV/c, and 
less than 1.0% above 2.3GeV/c. 

The systematic uncertainties inherent in the model- 
ing of the signal decays to non-resonant hadronic states 
and their impact on the signal yield have been studied 
by varying the SF parameters. We try four sets of SF 
parameters, two derived from the recent analysis of the 
B — > Xj7 decays [25] based on the semi-inclusive pho- 
ton spectrum and moments derived from this spectrum, 
one derived from moments in inclusive B — > XJtv de- 
cays [26] , and one from a combined fit [27] to all moments 
measured by the BABAR Collaboration. For each set of 
SF parameters we calculate the signal momentum spec- 
trum and repeat the fit to the data. We observe small 
changes in the fitted B — > X c lv background which re- 
sult in changes of the signal yield. Taking into account 
the errors and correlations of the measured SF param- 
eters, we derive the errors listed separately in Table II. 
The impact is largest for the signal regions extending to 
lower momenta, where this becomes the largest source of 
systematic error. 

Not included in this estimate is the sensitivity of the 
signal yield to the event selection criteria, specifically 
those based on the variables IZ2 and p m iss- These se- 
lection criteria influence not only the signal, but more 
so the background distributions. Details are discussed 
below. 



B. Non-BB Background 

Systematic errors in the subtraction of the non-BB 
background could be introduced by the choice of the fit- 
ting function describing the electron spectrum and by the 
uncertainty in the relative normalization of the on- and 
off-resonance data. 

To assess the uncertainty in the shape of this back- 
ground we have compared fits with different parameteri- 
zations of the fit function. In addition to the exponential 
function described above (Eq. 2), we have tried linear 
combinations of Chebyshev polynomials up to fifth or- 
der. The resulting fits are equally consistent with the 
data. The differences in the non-£?_B background esti- 
mates between different parameterizations are less than 
0.5%. Above 2.8 GeV/c the number of observed events 
in the on-resonance data sample agrees to 0.3% with 
the number of events predicted from the fit to the off- 
resonance sample. 

If the relative normalization is treated as a free fit pa- 
rameter, its deviation from the value based on luminosity 
measurements is less than one standard deviation. Thus, 
we use the more accurate value based on luminosity mea- 
surements. As a systematic error for the non-BB back- 
ground we take 0.5% of this background contribution, 
which includes the errors of the normalization factor and 
the background shape approximation. 
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C. B — > J/ipX Background 

J/ip decays to e + e~ pairs are vetoed by a restriction 
on the di-electron invariant mass. However, this veto is 
only about 50% efficient, primarily because of acceptance 
losses. The remaining, mostly single-electron background 
is estimated from simulation. We observe a difference 
of (5.0 ± 2.7)% between the veto efficiency for electron 
pairs in data and simulation, and thus assign a 5% error 
to the residual background. This background amounts 
to 18% of the signal for p e > 2.0 GeV/c and 10% for 
p e > 2.3 GeV/c and the resulting uncertainty on the sig- 
nal branching fraction is estimated to vary from 0.9% to 
0.5%. The background from ip{2S) — > e+e~ decays is 
significantly smaller, and thus its uncertainty is negligi- 
ble. 



D. BB Background 

The shapes of the BB backgrounds are derived from 
MC simulations. The branching fractions for exclusive 
scmilcptonic B — ► X c ev decays are currently not pre- 
cisely known. Thus the electron spectra from inclusive 
B — > X c ev decays may differ from those of the simula- 
tion. For this reason, we have introduced scale factors 
in the fits to the electron spectrum to adjust the relative 
normalization of the various contributions. To test the 
sensitivity to the shape of the dominant contributions, 
we have varied the form factors for decays to D*ev and 
Dev, and changed the relative proportion of contribu- 
tions from narrow and wide resonances to D**ev decays. 

For B — *■ Dev and B — *■ D*ev decays we use HQET 
parameterizations [36, 37] of the form factors. To study 
the impact of the uncertainties in the measured form fac- 
tors, we reweight the MC-simulated spectrum for a given 
decay mode to reproduce the change in the spectrum due 
to variations of the form-factor parameters, and repeat 
the standard fit to the data. From the observed changes 
in the signal yield as a function of the choice of the form- 
factor parameters for D*ev decays, we assess the system- 
atic error on the signal yield by taking into account the 
measured form- factor parameters, p 2 , R\, and R2, their 
errors, and their covariance matrix [40]. For Dev decays, 
we rely on a measurement of p 2 D by the CLEO [38] and 
Belle Collaborations [39]. Similarly, we estimate the im- 
pact of the uncertainty in p 2 D by comparing the default 
fit results with spectra corresponding to variations of p 2 D 
by one standard deviation. We take the shift of the signal 
yield as a systematic error. 

To assess the impact of the poorly known branching 
fractions for various D**ev decay modes on the shape 
of the electron spectrum, we have repeated the fits with 
the relative branching fractions for the individual decay 
modes changed by up to 50%. As long as we do not elim- 
inate the decays to the two narrow resonances, £>i(2437) 
and i?2(2459), we obtain consistent results. Specifically, 
if we eliminate the decays involving the two wider res- 



onances, Z?o(2308) and D[(2A60), the results change by 
less than 3%. We adopt this change as the estimate of 
the systematic error due to the uncertainty of decays to 
D** states. 

Similarly, we vary the branching fractions for sec- 
ondary electrons from semileptonic D decays by 10% and 
adopt the observed change as a systematic error. In ad- 
dition, there is a small contribution from events which 
contain a B — > X u ev decay, but contribute to the back- 
ground rather than the signal, because the track iden- 
tified as a signal electron does not originate from this 
decay. We estimate the uncertainty of this very small 
contribution to be 30%. 

For background from hadronic B decays, the uncer- 
tainty in the spectrum is primarily due to the uncer- 
tainty in the momentum-dependent hadron misidentffica- 
tion. The uncertainties of misidentification probabilities 
are estimated to be 20% and 30%, for pions and kaons, re- 
spectively. The uncertainty in the fractions of pions and 
kaons is taken as the difference between simulated and 
observed charged particle spectra, which is about 5% for 
pions and kaons. With these uncertainties in the hadron 
misidentification backgrounds, the fractional error in the 
number of subtracted background events is ~ 20% for 
pions and ~ 30% for kaons. In addition, there is a small 
background from protons and from unidentified particles; 
its total uncertainty is estimated to be about 50% smaller 
than for identified kaons. 



E. Uncertainty in the B Meson Momentum 
Spectrum 

The non-zero momentum of the B meson in the T(4S*) 
rest frame affects the shape of the electron spectrum near 
the endpoint. To estimate the systematic error associated 
with the uncertainty in the initial B-meson momentum 
spectrum, we compare the simulated and measured en- 
ergy spectra for fully reconstructed charged B mesons for 
different data taking periods. The widths of the energy 
distributions agree well for all data, but in some of the 
data sets we observed a shift in the central value of up 
to 2.2 MeV relative to the simulation, which assumes a 
fixed center-of-mass energy. We correct the simulation 
for the observed shifts, and assess the effect of the uncer- 
tainty of 0.13 MeV in this shift on the branching fraction 
measurement. 



F. Bremsstrahlung and Radiative Corrections 

For comparison with other experiments and with the- 
oretical calculations, the signal spectrum resulting from 
the fit is corrected for bremsstrahlung in the detector and 
for final-state radiation. Corrections for QED radiation 
in the decay process are simulated using PHOTOS [31]. 
This simulation includes multiple-photon emission from 
the electron, but does not include electroweak correc- 
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tions for quarks. The accuracy of this simulation has 
been compared to analytical calculations performed to 
0(a) [31]. Based on this comparison we assign an uncer- 
tainty of 20% to the PHOTOS correction, leading to an 
uncertainty in the signal yield of about 1%. 

The uncertainty in the energy loss of electrons due to 
bremsstrahlung in the beam pipe and tracking system 
is determined by the uncertainty in the thickness of the 
detector material, estimated to be (0.0450 ± 0.0014)X 
at normal incidence. The thickness of the material was 
verified using electrons from Bhabha scattering as a func- 
tion of the polar angle relative to the beam. The impact 
of the uncertainty in the energy loss on the signal rate 
was estimated by calculating the impact of an additional 
0.0014X of material. 



G. Sensitivity to the Event Selection 



X u ev according to 

2e{Ap)N BE 

Here Ntot refers to the total number of electron candi- 
dates selected in the on-resonance data and Nb g refers 
to the total background, from non-BB and BB events, 
as determined from the fit to the spectrum. e(Ap) is 
the total efficiency for selecting a signal electron from 
B — > X u ev decays (including bremsstrahlung in the de- 
tector material), and S ra d accounts for the distortion of 
the electron spectrum due to final-state radiation. This 
is a momentum-dependent correction, derived from the 
MC simulation based on PHOTOS [31]. The total num- 
ber of produced BB events is N BB - = (88.36 ± 0.02 s t at ± 
0.97 syst ) x 10 6 . 



We have checked the sensitivity of the fits to the elec- 
tron spectrum to changes in the event selection. We have 
also assessed the impact of the momentum-dependent un- 
certainty in the electron efficiency on the fitted signal 
yield. These variations of the event selection change the 
signal efficiency and lead to variations of up to 50% in 
the size of the non-BB background, and up to 20% in 
the BB background. 

Though some of the observed changes in the efficiency- 
corrected signal yield may already be covered by the 
form- factor and other variations, we conclude that these 
tests do reveal significant changes that have to be ac- 
counted for. 

The largest variation (5%) is observed for changes in 
the restriction on ratio of the Fox- Wolfram moments, 
IZ2, from the default value of 0.5 to 0.6. Other sizable 
variations are observed for changes in the restrictions on 
the absolute value and direction of the missing momen- 
tum vector. IZ2 and the missing momentum are quan- 
tities that are derived from the measured momenta of 
all charged and neutral particles in the event, and are 
therefore sensitive to even small differences in data and 
simulation. We interpret the observed changes as repre- 
sentative for the uncertainties in the MC simulation of 
the selection of signal and background events and adopt 
the observed changes between the default fits and the fits 
with looser selection criteria as systematic errors. Adding 
the observed changes in quadrature leads to a relative 
systematic error of between 5% and 8% on the partial 
branching fraction. 



V. RESULTS 

Determination of the Partial B 
Branching Fraction 



X u eu 



> 
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For a given interval Ap in the electron momentum, we 
calculate the inclusive partial branching fraction B — > 



FIG. 5: The differential branching fraction for charmless 
semileptonic B decays (data points) as a function of the 
electron momentum (in the T(4S r ) rest frame) after back- 
ground subtraction and corrections for bremsstrahlung and 
final state radiation, compared to the Monte Carlo simula- 
tion (histogram). The errors indicate the statistical errors on 
the background subtraction, including the uncertainties of the 
fit parameters. For the signal simulation, the SF parameters 
are extracted from a combined fit to all BABAR moments. 



The differential branching fraction as a function of the 
electron momentum in the T(4S I ) rest frame is shown 
in Fig. 5, fully corrected for efficiencies and radiative ef- 
fects. The data are well reproduced by the signal simu- 
lations using the SF parameters derived from the com- 
bined fit to all moments measured by the BABAR Collab- 
oration [27], specifically mf F (1.5GeV) = 4.59GeV/c 2 , 
/4 SF (1.5GeV) = 0.21 GcV 2 . The partial branching frac- 
tions for the five overlapping electron momentum inter- 
vals are summarized in Table III. The stated errors on 
AB represent the statistical and total systematic uncer- 
tainties of the measurement, including the uncertainty 
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due to the sensitivity to the SF parameters, as stated 
in Table II. As the lower limit on the electron momen- 
tum decreases, the statistical and systematic errors are 
more and more dominated by the B — > X c ev background 
subtraction. 



B. Extraction of the Total Charmless Branching 
Fraction and |Vub| 

As mentioned earlier, we use two sets of theoretical cal- 
culations to extract \V u b\ from the partial electron spec- 
trum. The first, and so far the most commonly used, 
method derives \V u b\ from the total charmless scmilep- 
tonic branching fraction and the average B lifetime, 
n = (1.604 ± 0.012) ps [43], as follows, 

X (1.0 =b 0.028 pert+nonpert 

±0.039 m J. (5) 

Here the first error represents the linear sum of the uncer- 
tainties of the perturbative and non-perturbative QCD 
corrections, and the second error is due to the uncer- 
tainty in mi,. An overall correction of 0.7% is included 
to account for QED corrections. This formulation [44- 
46] has been updated to take into account the recent 
measurement [26] of nib, and other parameters of 
the heavy quark expansion in the kinetic mass scheme, 
specifically m£ m (1.0GeV) = (4.61 ± 0.07) GeV/c 2 and 

/4 (fcin) (1.0GeV) - (0.45 ± 0.05) GeV 2 . 

We determine the total branching fraction, 

B(B -> X u ev) = AB(Ap)/f u (Ap), (6) 

where f u (Ap) is the fraction of the electron spectrum in 
a given momentum interval Ap. The values of f u (Ap) 
are estimated based on the DN calculations [17] using 
the exponential parameterization of the SF, with the 
SF parameters extracted from fits to the photon spec- 
trum in semi-inclusive B — > X s ~f decays, as measured by 
the BABAR Collaboration [25], A SF = (0.49±g;Jg) GeV/c 2 , 
\f F = (-0.24+Q °g) GeV 2 , with a correlation coefficient 
of —0.94. We obtain very similar results for the two other 
functional forms suggested to describe the SF [17]. 

The results for the predicted fraction f u (Ap), the total 
charmless branching fraction, B, and |V^&| are presented 
in Table III. The first error on /„ refers to the experi- 
mental error of the SF parameters from the measurement 
of the inclusive photon spectrum. It includes the uncer- 
tainty of the background subtraction and the extrapola- 
tion to decays to unmeasured X s states. We have taken 
into account the stated error of the SF parameters, in- 
cluding their correlation. Specifically, we have taken as 
an error on f u (Ap) the maximum deviation of the f u (Ap) 
from its central value for selected values of the SF param- 
eters on the error ellipse. The second error accounts for 
the dependence on the a s scale, for the uncertainty in 



the form of the SF, and for the uncertainty in the theo- 
retical prediction of /„ from the B — > X s j measurement. 
As suggested by M. Neubert [47] this error has been es- 
timated by varying the values of A SF and Xf F by 10%. 

The errors listed for B and \ V u b \ are specified as follows. 
The first error reflects the error on the measurement of 
AB, which includes statistical and experimental system- 
atic uncertainties, except for the uncertainty in the SF 
parameters. The second error is due to experimental un- 
certainty of SF parameters affecting both f u (Ap) and 
AB. The third error is the theoretical uncertainty of 
f u (Ap). The fourth error on |V^&| accounts for the the- 
oretical uncertainty in the translation from B to |V^f,|, 
as specified in Eq. 5. This error also depends on the b- 
quark mass and thus is correlated with the theoretical 
uncertainty on the SF. 

The results for the total branching fraction B and \ V u b\ 
obtained from the different momentum intervals are con- 
sistent within the experimental and theoretical uncer- 
tainties. For intervals extending below 2.3GeV/c, the to- 
tal errors on B and | V u b | do not depend very strongly on 
the chosen momentum interval. While the errors on AB 
are smallest above the kinematic endpoint for B — ► X c ev 
decays, the dominant uncertainty arises from the deter- 
mination of the fraction /„ and increases substantially 
with higher momentum cut-offs. The stated theoretical 
errors on /„ , acknowledged as being underestimated [47] , 
do not include uncertainties from weak annihilation and 
other power-suppressed corrections. Assuming that one 
can combine the experimental and theoretical errors in 
quadrature, the best measurement of the total branch- 
ing fraction is obtained for the momentum interval 2.0- 
2.6 GeV/c. 

Though the BABAR measurement of the photon spec- 
trum [25] results in the best estimate for the SF param- 
eters, we have also considered sets of SF parameters ob- 
tained from photon spectra measured by the CLEO [48] 
and Belle [49] Collaborations. These parameters are 
listed in Table IV. In Table V the results obtained 
for these different SF parameters based on the BABAR 
scmilcptonic data and on the DN calculations are listed 
for the momentum interval 2.0 — 2.6 GeV/c. The differ- 
ences between the SF parameters obtained by the CLEO 
and Belle Collaborations and the BABAR results are com- 
parable to the experimental errors on these parameters. 
These differences affect the signal spectrum, and thereby 
the fitted background yield. The effect is small for high 
momentum region and increases for the signal intervals 
extending to lower momenta. The impact of the SF pa- 
rameters on the partial branching fractions is included in 
the total error (see Table II). 

The second method for extracting \V u b\ is based on 
recent BLNP calculations [24]. In this framework the 
partial branching fraction AB is related directly to \V u b\: 

where C(Ap) is the prediction for the partial rate for B — > 
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TABLE III: The partial (AS) and total (B) branching fraction for inclusive B — > X u ev decays and \V u b\ for five electron 
momentum intervals. The spectral fractions f u are determined using SF parameters extracted from a fit to the photon 
spectrum in B — » X s y decays [25] based on calculations by DeFazio and Neubert [17] and Kagan and Neubert [18]. The errors 
are explained in the text. 

Ap(GeV/c) Ag(HT 3 ) f u (Ap) 6(l(r 3 ) |K o |(10~ 3 ) 

2.0 - 2.6 0.479 ± 0.033 ± 0.050 0.298 ± 0.029 ± 0.015 1.61 ± 0.18 ± 0.25 ± 0.08 3.80 ± 0.21 ± 0.29 ± 0.10 ± 0.18 

2.1 - 2.6 0.350 ± 0.020 ±0.033 0.222 ± 0.026 ± 0.016 1.58 ± 0.16 ± 0.24 ± 0.11 3.77 ± 0.19 ± 0.28 ± 0.13 ± 0.18 

2.2 - 2.6 0.231 ±0.010 ±0.018 0.149 ± 0.020 ± 0.016 1.55 ± 0.14 ± 0.23 ± 0.16 3.73 ± 0.17 ± 0.27 ± 0.19 ± 0.18 

2.3 - 2.6 0.146 ±0.006 ±0.010 0.086 ± 0.013 ± 0.013 1.71 ± 0.13 ± 0.25 ± 0.26 3.92 ± 0.15 ± 0.29 ± 0.30 ± 0.19 

2.4 - 2.6 0.075 ±0.004 ±0.006 0.039 ± 0.006 ± 0.009 1.95 ± 0.20 ± 0.32 ± 0.45 4.18 ± 0.21 ± 0.35 ± 0.48 ± 0.20 



TABLE IV: SF parameters (at a scale of 1.5 GeV) measured 
by different experiments, based on two different theoretical 
calculations, top: DN [17, 18], bottom: BLNP [24]. 



Experiment SF Input A (GeV/c 2 ) Ai(GeV 2 ) 



BABAR 


(spectrum) X s -y 


[25] 


0.491°;^ 


-0.24in 9 8 


CLEO 


(spectrum) X s ~{ 


[48] 


0.54«; 2 * 


-0.34t°il 


Belle 


(spectrum) Xs') 


[49] 


0.66i° ;» 




Experiment SF Input A (GeV/c 2 ) 


A ( GeV 2 ) 


BABAR 


(spectrum) X s j 


[25] 


n 6i + 07 

u - D1 -0.07 


o.i6lg:SS 


BABAR 


(moments) X a ~{ 


[25] 


U - '°-0.13 


o.35±g : iS 


BABAR 


(moments) X c iv 


[26] 


0.67 ±0.08 


0.15 ±0.07 


BABAR 


(combl. moments) 


[27] 


0.69 ±0.05 


0.21 ±0.05 



X u ev decays (in units of ps _1 ). In these calculations the 
leading order SF is constrained by the HQE parameters, 
obtained either from the B — > X s j or B — > X c ev decays, 
or both. 

The values of the SF parameters extracted from the 
BABAR analyses of inclusive B — -> X s ^ [25], B — > 
X c ev [26] decays, and the combined fit [27] to all mo- 
ments measured by the BABAR Collaboration are listed 
in Table IV. Note that the definitions of shape func- 
tions and the SF parameters are different for the DN 
and BLNP calculations. The different SF parameters 
and their measurement errors are also shown in Fig. 6. 

The SF parameters based on B — > X s ~f data only 
are extracted from either a fit to the photon spectrum 
or to the first and second moments of this spectrum 
in the "shape function" scheme. The HQE parame- 
ters extracted from fits to measured moments in the ki- 
netic mass scheme have been translated into the "shape 
function" scheme at the appropriate scale. Specifically, 
the HQE parameters extracted from the moments in 
B — > X c ev decays have been translated based on two- 
loop calculations [23]. The HQE parameters resulting 
from the combined fit to moments of the photon, lepton, 
and hadron mass spectra in the kinetic scheme are used to 
predict the first and second moments of the photon spec- 
trum down to photon energies of 1.6 GeV, based on cal- 
culations by Benson, Bigi, and Uraltsev [50]. The lower 
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FIG. 6: (color online) The fitted values and contours corre- 
sponding to A^ 2 = 1 for the four sets of SF parameters (see 
Table IV) based on the calculations of BLNP, extracted from 
the photon energy spectrum (short dash - red) and from the 
photon energy moments (dot-dash - green) in B — » Xs'y, from 
the lepton energy and hadron mass moments in B — » X c ev 
decays (long dash - black), as well as from the combined fit 
to moments (solid - blue) measured by the BABAR Collabo- 
ration. Also shown are two straight lines indicting values of 
the SF parameters, for which the partial branching fraction 
(dotted - magenta) and |V^b| (solid - light blue) are constant. 



limit on the photon energy is chosen such that the es- 
timated cut-induced perturbative and non-perturbative 
corrections to the HQE are negligible. From these pre- 
dicted moments, the SF parameters are extracted us- 
ing the next-to-leading order calculations in a framework 
that is consistent with the one used for the determination 
of | Kb | [24]. 

The smallest errors on the SF parameters are obtained 
from the fit to the photon spectrum and the combined fit 
to all moments. The fit to the photon spectrum is most 
sensitive to the high end of the photon energy spectrum, 
and relies on the theoretical prediction for the shape of 
the spectrum down to low photon energies. Since this 
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TABLE V: Comparison of measurements of the partial (AB) and total (B) branching fraction for inclusive B — > X u ev decays 
and \ V u b \ for the electron momentum interval 2.0 to 2.6 GeV/c. The results are obtained for SF parameters (listed in Table IV) 
extracted from different experiments. The first three measurements are based on DN [17, 18] calculations, the remaining four 
on BLNP [24] calculations, based on SF parameters extracted from the photon spectra and energy moments, the B — ► X c lv 
moments [26], and a combined fit to moments [27]. The errors are explained in the text. 



Experiment 


SF input 


AB(1CT 3 ) 


23(10~ 3 ) 


IKbKio -3 ) 


BABAR 


Xs'j (spectrum) 


0.479 ± 0.033 ± 0.050 


1.61 ±0.18 ±0.25 ±0.08 


3.80 ± 0.21 ± 0.29 ± 0.10 ± 0.18 


CLEO 


Xs'y (spectrum) 


0.491 ±0.036 ±0.061 


1.75 ±0.20 ±0.48 ±0.11 


3.97 ± 0.23 ± 0.54 ± 0.12 ± 0.19 


Belle 


X s 7 (spectrum) 


0.548 ± 0.038 ± 0.057 


2. 24 ±0.25 ±0.27 ±0.20 


4.48 ± 0.25 ± 0.27 ± 0.20 ± 0.22 


BABAR 


Xa'j (spectrum) 


0.514 ± 0.037 ± 0.055 


1.81 ±0.20 io.M ±0.11 


3.80 ±0.21 toil ±°- 18 


BABAR 


Xs'j (moments) 


0.577 ±0.041 ±0.082 


2.57 ±0.29 lo.gg ±0.23 


4.86 ± 0.28 tlil ± 0.26 


BABAR 


X c iu (moments) 


0.569 ± 0.039 ± 0.090 


2.17 ±0.24 ±0.15 


4.30 ± 0.24 toll ± °- 21 


BABAR 


combined fit to moments 


0.572 ± 0.041 ± 0.065 


2.27 ±0.26 to.fs ±0.17 


4.44 ± 0.25 to ll ± 0.22 



shape is not directly calculable, several forms of the SF 
are used to assess the uncertainty of this approach. The 
use of two sets of the first and second moments of the 
photon spectrum, above 1.90 and above 2.09 GeV, is less 
powerful, due to much larger statistical and systematic 
errors, but insensitive to the theoretical knowledge of the 
detailed shape of the spectrum. The SF parameters ob- 
tained from moments of the photon spectrum above 1.90 
GeV/c agree with those obtained from the global fit to 
the moments, but also have larger errors. Nevertheless, 
the inclusion of the photon energy moments significantly 
improves the sensitivity of the global fit to more than 30 
measured moments. 

The results for the partial branching fractions AB and 
| V u b | based on the BLNP calculations are listed in Ta- 
bles VI, VII, VIII, and IX for the four sets of SF param- 
eters. 

The errors cited in these tables are denned and deter- 
mined in analogy to those in Table III. The first error on 
the predicted rate £ accounts for the uncertainty due to 
the errors in measured parameters of the leading SF, the 
second error refers to the theoretical uncertainties in the 
subleading SFs, and variations of scale matching, as well 
as weak annihilation effects. For \V u b\, the first error is 
the experimental error on the partial branching fraction, 
which includes the statistical and the experimental sys- 
tematic uncertainty, the second error includes systematic 
uncertainties on the partial branching fraction and ( due 
to the uncertainty of the SF parameters, and the third 
error is the theoretical uncertainty on £, estimated using 
the prescription suggested by BLNP. 

In Table V the results obtained for these different SF 
parameters based on the BABAR semileptonic data and 
on the BLNP (and DN) calculations are listed for the 
momentum interval 2.0 — 2.6 GeV/c. The observed differ- 
ences are consistent with the total error stated; they are 
largest for the SF parameters extracted from the fit to 
the photon spectrum as compared to the moments of the 
photon spectrum. 

For all four sets of SFs we observe a tendency for the 



total branching fraction, and therefore also \V u b\, to be 
slightly larger at the higher momentum intervals, but the 
uncertainties in the predicted rates £ are very large for 
the highest momentum interval. 

Based on the BLNP calculations [24] of the inclusive 
lepton spectra, we have also determined the total B — > 
X u ev branching fraction. The results are presented in 
Table X. 

The results for \V u b\ extracted for the BLNP calcula- 
tions are close to those obtained for the DN calculations 
(see Table V). In fact, the results based on the fit to 
the photon spectrum measured by the BABAR Collabo- 
ration are identical for all electron momentum ranges, 
even though the partial branching fractions differ by one 
standard deviation of the experimental error (sec Tables 
III and VI). Changing the ansatz for the SF from the 
exponential to a hyperbolic function [24] has no impact 
on the results. 



VI. CONCLUSIONS 

In summary, we have measured the inclusive electron 
spectrum in charmless semileptonic B decays and derived 
partial branching fractions in five overlapping electron 
momentum intervals close to the kinematic endpoint. 
We have extracted the partial and total branching frac- 
tions and the magnitude of the CKM element \V u b\ based 
on two sets of calculations: the earlier ones by DeFazio 
and Neubert [17] and Kagan and Neubert [18], and the 
more comprehensive calculations by Lange, Neubert and 
Paz [24], as summarized in Table V. Within the stated 
errors, the measurements in the different momentum in- 
tervals are consistent for both sets of calculations. 

We adopt the results based on the more recent calcula- 
tions (BLNP) [24], since they represent a more complete 
theoretical analysis of the full electron spectrum and re- 
late the SF to the HQE parameters extracted from inclu- 
sive B — > X s j and B — > X c lv decays. We choose the SF 
parameters obtained from the combined fit to moments 
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TABLE VI: The partial branching fraction AB, the predicted partial rate £ for B — » X u eu decays and \V u b\ for five 
electron momentum intervals, using the SF parameters from the photon spectrum in semi-inclusive B — * X 3 *y decays, 
(A SF = (0.611^7) GeV/c 2 , ^ 2SF = (0.16tg;£g) GeV 2 , fit to spectrum) [25] based on BLNP calculations [24]. The errors 
are explained in the text. 

Ap(GeV/c) AB (1CT 3 ) ((Ap)^ 1 ) |K b |(10~ 3 ) 

2.0-2.6 0.514 ±0.037 ±0.055 22.2 ±3.9 ±2.1 3.80 ± 0.21 t°H ± 0.18 

2.1 - 2.6 0.366 ±0.021 ±0.034 16.1 ±3.3 ± 1.7 3.76 ± 0.20 tg;fjj ± 0.20 

2.2 - 2.6 0.236 ±0.011 ±0.019 10.5 ±2.4 ± 1.5 3.75 ± 0.17 ±g;^ ± 0.27 

2.3 - 2.6 0.147 ±0.006 ±0.010 5.8 ±1.5 ± 1.6 3.98 ± 0.16 ±g;|| ± 0.55 

2.4 - 2.6 0.075 ± 0.004 ± 0.006 2.5 ±0.7 ±2.1 4.32 ± 0.22 ±° H ± 1.81 



TABLE VII: The partial branching fraction AB, the predicted partial rate £ for B — > X u ev decays and \V u b\ for five electron 
momentum intervals, using the SF parameters from the moments of the photon energy spectrum in semi-inclusive B — > X s j 
decays, (A SF = (0.75±%\l) GeV/c 2 , ^ 2SF = (0.35±q\\\) GeV 2 , fit to moments) [25] based on BLNP calculations [24]. The 
errors are explained in the text. 

Ap(GcV/c) AB (10~ 3 ) C(A P )(ps- 1 ) |14b|(10~ 3 ) 

2.0 - 2.6 0.577 ± 0.041 ± 0.082 15.2 ±5.1 ± 1.6 4.86 ± 0.28 t, 1 ,;^ ± 0.26 

2.1 - 2.6 0.392 ±0.024 ±0.043 10.5 ±4.1 ± 1.5 4.82 ± 0.25 ±J |g ± 0.33 

2.2 - 2.6 0.243 ±0.011 ±0.021 6.5 ± 3.0 ± 1.4 4.82 ± 0.22 ±J || ± 0.52 

2.3 - 2.6 0.148 ±0.006 ±0.010 3.5 ±1.8 ±1.6 5.15 ± 0.20 ± 1.20 

2.4 - 2.6 0.075 ±0.004 ±0.006 1.5 ±0.8 I 2 ; 3 5.62 ± 0.29 t\f 4 ± 4.27 



of inclusive distributions measured by the BABAR Collab- 
oration rather than the single most precise measurement 
of the SF parameters obtained from the recent BABAR 
measurement [25] of the semi-inclusive photon spectrum 
in B — > X s ^f decays. Assuming it is valid to combine 
the experimental and the estimated theoretical errors in 
quadrature, and taking into account the fraction of the 
signal contained in this interval, we conclude that the 
best measurement can be extracted from the largest mo- 
mentum interval, 2.0 to 2.6GeV/c. For this momentum 
interval the partial branching fraction is 

AB(B -» X u ev) = (8) 
(0.572 ± 0.041 stat ± 0.065 syst ) x lO" 3 . 

Here the first error is statistical and the second is the to- 
tal systematic error, as listed in Table II. In addition to 
the systematic uncertainty due to the signal extraction, 
the normalization, and various small corrections, this 
error also includes the observed dependence of the ex- 
tracted signal on the choice of the SF parameters. Based 
on the BLNP method, we obtain a total branching frac- 
tion of 

B(B -> X u ev) = (9) 
(2.27±0.26 e:cp t° fs sf ±0.17 theory )xW- 3 , 

and 

\V ub \ = (10) 
(4.44±0.25 exp ±gil sf ± 0-22 theory ) x lO" 3 . 



Here the first error represents the total experimental un- 
certainty, the second refers to the uncertainty in the SF 
parameters from the combined fit to moments, and the 
third combines the stated theoretical uncertainties in the 
extraction of \V u b\, including uncertainties from the sub- 
leading SFs, weak annihilation effects, and various scale- 
matching uncertainties. No additional uncertainty due to 
the theoretical assumption of quark-hadron duality has 
been assigned. 

The improvement in precision compared to earlier 
analyses of the lepton spectrum near the kinematic end- 
point can be attributed to improvements in experimen- 
tal techniques, to higher statistics, and in particular, to 
improved background estimates, as well as significant ad- 
vances in the theoretical understanding of the SFs and ex- 
traction of the SF parameters from inclusive spectra and 
moments. While earlier measurements were restricted 
to lepton energies close to the kinematic endpoint for 
B — > X c lv decays at 2.3 GeV/c and covered only 10% 
of the B — > X u lv spectrum, these and other more recent 
measurements have been extended to lower momenta, in- 
cluding about 25% of the spectrum, and thus have re- 
sulted in a significant reduction in the theoretical uncer- 
tainties on \y u b\. 

The determination of \V u b\ is currently limited primar- 
ily by our knowledge of SF parameters. An approximate 
linear dependence of \V u b\ on these parameters is 

^ = 1.31^-0.04^. (11 ) 
\Vub\ A a4 
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TABLE VIII: The partial branching fraction AB, the predicted partial rate £ for B — » X u eu decays and \V u b\ for five electron 
momentum intervals, based on the SF parameters from hadron mass and lepton moments in B — > X c ev decays (A SF = 
(0.67 ± 0.08) GeV/c 2 , ni SF = (0.15 ± 0.07) GeV 2 ) [26] based on BLNP calculations [24]. The errors are explained in the text. 



Ap ( GeV/c) 



AB(10" 



C(Ap)(p S " 1 ) 



|Kf,|(10- 



2.0 - 2.6 

2.1 - 2.6 

2.2 - 2.6 

2.3 - 2.6 

2.4 - 2.6 



0.569 ±0.039 ±0.090 
0.391 ± 0.022 ± 0.048 
0.243 ±0.011 ±0.022 
0.148 ±0.006 ±0.010 
0.075 ± 0.004 ± 0.006 



19.2 ±3.9 ±1.9 
13.5 ±3.3 ±1.6 
8.3 ± 2.6 ± 1.5 
4.3 ± 1.7 ± 1.7 



1.7 ±0.8 



+2.4 



4.30 ±0.24 toll ±°- 21 

4.26 ±0.22 toll ±0.25 

4.27 ±0.19 toll ±0-37 
4.65 ±0.18 to 9 7 l ±0.91 

-1.29 
-1.00 



TABLE IX: The partial branching fraction AB, the predicted partial rate £ for B — > X u ev decays and \V u b\ for five electron 
momentum intervals, based on the SF parameters from the combined fit to BABAR moments (A = (0.69 ± 0.05) GeV/c 2 , 
H 2SF = (0.21 ± 0.05) GeV 2 ) [27] based on BLNP calculations [24]. The errors are explained in the text. 



Ap ( GeV/c) 



AB(l(r 



C(Ap)(ps- 1 ) 



\v ub \(io- 



2.0 - 2.6 

2.1 - 2.6 

2.2 - 2.6 

2.3 - 2.6 

2.4 - 2.6 



0.572 ± 0.041 ± 0.065 
0.392 ±0.023 ±0.038 
0.243 ±0.011 ±0.020 
0.148 ±0.006 ±0.010 
0.075 ± 0.004 ± 0.006 



18.1 ±2.3 ±1.8 
12.6 ±1.9 ±1.5 
7.8 ± 1.4 ± 1.4 
4.1 ±0.9 ± 1.7 
1.7 ±0.5 tf.7 



4.44 ± 0.25 t H ± 0.22 
4.40 ± 0.23 to.tl ± 0- 27 
4.40 ± 0.20 toil ± 0-41 
4.74 ±0.18 t°H ±0.96 
5.29 ±0.27 tilt ±3.66 



TABLE X: The total (B) branching fraction for inclusive 
B — > X u ev decays for five electron momentum intervals. The 
spectral fractions f u are based on calculations by Lange, Neu- 
bert and Paz [24] using SF parameters extracted from the 
combined fit [27] to all BABAR moments. The error defini- 
tions are the same as in Table III, and they are explained in 
the text above. 



Ap ( GeV/c) 



MM 



2.0 - 2.6 0.252 ± 0.018 ± 0.019 2.27 ± 0.26 ±g-_|| ± 0.17 

2.1 - 2.6 0.176 ± 0.017 ± 0.019 2.22 ± 0.23 t°ir ± 0.24 

2.2 - 2.6 0.109 ± 0.014 ± 0.020 2.22 ± 0.20 t°H ± 0.40 

2.3 - 2.6 0.057 ± 0.009 ± 0.023 2.58 ± 0.20 t°f 6 ± 1-04 

2.4 - 2.6 0.023 ± 0.005 ± 0.033 3.21 ± 0.33 t° H tUl 



for A = 0.69 GeV/c 2 and /4 = 0.21 GeV 2 . Thus the un- 
certainty on the 6-quark mass dominates. It should be 
noted that this dependence on A is a factor of two smaller 
for measurements based on the DN calculations. 

These results are in excellent agreement with earlier 
measurements of the inclusive lepton spectrum at the 
T(45 f ) resonance, but their overall precision surpasses 
them [8-10, 12]. The earlier results were based on the 
DN calculations. We observe that for the same experi- 
mental input, i.e. the same measured lepton and photon 
spectra, the extracted values of \V u b\ based on DN cal- 
culations agree very well with those based on BNLP cal- 
culations for the various momentum ranges under study, 
even though the corresponding partial branching frac- 



tions may differ by one standard deviation. 

The results presented here are also comparable in 
precision to, and fully compatible with, inclusive mea- 
surements recently published by the BABAR [51, 52] 
and Belle [53] Collaborations, based on two-dimensional 
distributions of lepton energy, the momentum transfer 
squared and the hadronic mass, with SF parameters ex- 
tracted from B — ► A"„7 and B — ► X c lv decays. 
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